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J0316+4328: a Probable "Asymmetric Double" Lens 
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ABSTRACT 

We report a probable gravitational lens J0316+4328, one of 19 candidate asymmetric 
double lenses (2 images at a high flux density ratio) from CLASS. Observations with 
the Very Large Array (VLA), MERLIN and the Very Long Baseline Array (VLBA) 
imply that J0316+4328 is a lens with high confidence, ft has 2 images separated by 
0'.'40, with 6 GHz flux densities of 62 mJy and 3.2 mJy. The flux density ratio of 
~ 19 (constant over the frequency range 6-22 GHz) is the largest for any 2 image 
gravitational lens. High resolution optical imaging and deeper VLBI maps should 
confirm the lensing interpretation and provide inputs to detailed lens models. The 
unique configuration will give strong constraints on the lens galaxy's mass profile. 

Key words: gravitational lensing - cosmology:miscellaneous. 



1 INTRODUCTION 

Galaxy-scale gravitational lenses are valuable tools for the 
study of galactic structure. The locations and brightnesses 
of the images provide constraints on the overall lens galaxy 
profile, while remaining anomalies and time variations 
in the brightnesses due to micro-lensin g probe substruc- 
ture o n scales from dark matter clumps ( Mao fc Schneider] 
1 19981; iDalal fc Kochanekl 120021 ) and satellite galaxies 
llMcKean et alj 120071 ') to stars |Schechter fc Wambsgansa 
2002). Lens images are sensitive to the surface densities at 
the locations where they form, generally at galactocentric 
radii of a few kpc for bright lens images. Measuring the over- 
all lens galaxy profile over a range of scales requires more 
information, which can be provided by stellar vel ocity dis- 
persions at radii < 1 kpc llKoopmans et al.l|2 006) or weak 
lensing at radii of tens of kpc ( Gavazzi et al.ll2007l ). Lenses 
with central images measure the density profile on many 
scales using strong lensing alone, giving strong constraints 
at radii of ~ 100 pc, including any contribution from a cen- 
tral super-massive black hole l|winn et all 120041 '). 

For lensed quasars, the best measurement of the lens 
galaxy profile can be done if the quasar is radio-loud. The 
size of the e mitting region is larger in the r adio than in opti- 
cal or X-ray l|Koopmans fc de Bruvnll2000l ). making the im- 
age brightnesses less susceptible to microlensing, while cen- 
tral images are strongly affected by dust absorption and con- 
fusion in the opti cal regime but no t in the radio frequency 
range 5-40 GHz (|Winn et al.l 120041 ). The large st sample of 
radio-loud lenses has been provided by CLASS (|Mvers et al.l 
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120031 ; iBrowne et all 120031 ). which found 22 instances of a 
bright radio source lensed by a foreground galaxy. 

A crucial step in CLASS was the rejection of double 
sources with a flux density ratio > 10 : 1 l|Mvers et al.ll20o"3l ) 
after the initial 8.4 GHz VLA snapshot. This was done for 
two reasons; for statistical completeness of the lens sample 
it was necessary to be sure that all secondary components 
could be detected reliably and it also significantly reduced 
the amount of follow-up observations required. A secondary 
compact component of similar brightness to the primary is 
likely to be a lensed image, while a secondary component 
much fainter than the primary is most likely to be weak 
structure belonging to the primary source. This strategy 
necessarily excluded "asymmetric doubles": 2 image lenses 
with a high flux density ratio. Such lenses are of interest 
beca use they are mos t likely to host observ able central im- 
ages (|Mao et al.ll200ll ; iBowman et al.ll2004l ). while even the 
absence of a central image gives the strongest constraints 
in an asymmetric doub le (|Rusin fc Mai l200ll ; iBovce et all 
2006; Zha ng et al.1 l2007h . Though the recognition of weak 
secondaries with ratios > 10 : 1 in CLASS may not have 
been 100% reliable, many are detectable. We have, therefore, 
conducted a program to identify such lenses by following up 
promising candidates which were initially rejected only on 
the basis of a high flux density ratio. 



2 SURVEY AND OBSERVATIONS 

We selected 18 candidates which were fitted by multiple 
compact components in the initial CLASS VLA snapshots, 
and for which the second brightest component was 10 — 30 
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Table 1. Summary of VLA observations. Flatter and steeper refer to the spectral indices of weaker components over the observed 
frequencies, relative to the primary compact source, e.g. J0008+6837 had 2 compact radio components, the fainter of these had a steeper 
radio spectrum between 5 and 8.5 GHz. Most of our candidates were typical radio sources with a flat-spectrum, compact core and 
a steeper spectrum, often extended, lobe or knot in a jet. Only J0316+3350, J0935+0719 and J0958+2948 could be considered lens 
candidates from the VLA observations. 



Source 




R.A. 




Dec. 


Observing Time 


(mins.) 


1\T ^4- ^r-i 

IN otcs 








(J2000) 






5 GHz 


8.5 GHz 


OO f ' IT,. 




J0008+6837 


00 


fW 
U o 


33.5 


68 


O / 


99 


11 


"1 1 
11 




Rejected: 2nd component steeper 


J0152+3350 


01 


52 


34.6 


33 


50 


33 


6 


8 




Rejected: 2nd component steeper 


J0316+3350 


03 


16 


50.9 


43 


28 


19 


11 


12 


17 


Possible Lens from VLA observations 


J0644+5955 


06 


41 


04.7 


59 


55 


21 


6 


7 




Rejected: 2nd and 3rd components extended, steeper 


J0812+4041 


08 


12 


03.0 


10 


11 


08 


18 


18 




Rejected: 2nd and 3rd components extended, steeper 


J0852+5922 


08 


52 


30.0 


59 


22 


50 


18 


19 




Rejected: 2nd component extended, steeper 


J0903+4651 


09 


03 


04.0 


46 


51 


04 


4 


5 


45 


Rejected: 2nd component extended 


J0935+0719 


09 


35 


01.1 


07 


19 


19 


12 


11 


23 


Possible Lens from VLA observations 


J0938+3934 


09 


38 


39.2 


39 


34 


20 


18 


18 




Rejected: 2nd component extended, steeper 


J0958+2948 


09 


58 


58.9 


29 


48 


04 


12 


13 


21 


Possible Lens from VLA observations 


J1131+5146 


11 


31 


16.5 


51 


46 


34 


18 


19 




Rejected: 2nd component extended, steeper 


J1252+1910 


12 


52 


27.8 


19 


10 


38 


18 


19 




Rejected: 2nd component extended, steeper 


J1253+6304 


12 


53 


17.5 


63 


04 


36 


18 


18 


24 


Rejected: 2nd component flatter, 3rd component extended 


J1343+2844 


13 


43 


00.2 


28 


44 


07 


6 


6 




Rejected: 2nd component extended, steeper 


J1540+1447 


15 


40 


49.5 


14 


47 


46 


4 


5 


28 


Rejected: 2nd component extended, steeper 


J1641+5115 


16 


41 


55.7 


51 


15 


47 


9 


10 




Rejected: 2nd component extended, steeper 


J2139+1027 


21 


39 


42.6 


10 


27 


43 


7 


7 


11 


Rejected: 2nd component steeper 


J2353+3231 


23 


53 


20.9 


32 


31 


44 


14 


14 




Rejected: 2nd component extended, steeper 



times fainter than the brightest. New observations of these 
sources were made with the VLA in A configuration at 5 
and 8.5 GHz on 2006 February 12 and 2006 February 19. 
Each source was observed for a few minutes in several blocks 
at different hour angles, giving better u-v coverage than the 
initial snapshots. Candidates with multiple compact compo- 
nents at 5 and 8.5 GHz were re-observed at 22 GHz on 2006 
March 20 or 2006 April 03. Most candidates were found to 
have an extended secondary component, or differing spectral 
indices between components, leaving 3 candidates as possi- 
ble lenses, J0316+4328, J0935+0719 and J0958+2948. The 
observations are summarised in Table Q] 

J0958+2958 has 2 compact radio sources separated by 
8", but these objects are distinct quasars with redshifts of 
2.06 4 and 2.744, rulin g out this system as a gravitational 
lens (|Lehar et al.ll200ll ). J0935+0719 has a series of radio ob- 
servations indicating that it might be a gravitational lens in 
which the source is a compact symmet ric object, and ismen- 
tioned as the candidate J0935+073 in iBrowne et all (|2003T ). 
We find that the flux density ratio is 19:1 at 5, 8.5 and 
22 GHz, supporting the lensing interpretation. A definitive 
answer on the nature of J0935+0719 requires better opti- 
cal data, both imaging and spectroscopic, and we do not 
consider this system further in this paper. 

This left J0316+4328 as a possible candidate lack- 
ing high resolution radio imaging. After completing the 
longer VLA observations we realised that another source, 
J0856+4935, looked promising in the initial CLASS snap- 
shot. On December 08 2006 and December 09 2006 we 
observed both sources at 6.0 GHz with the Multiple El- 
ement Radio-Linked Interferometric Network (MERLIN), 
spending 12 hours on each source. J0856+4935 had a dif- 
fuse secondary component and was rejected as a lens, while 
J0316+4328 showed 2 compact components and remained 
as a lens candidate with increased confidence. 
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Figure 1. 8.5 GHz VLA A array map of J0316+4328, made with 
super-uniform weighting. Components A and B are probably 2 
images of a lensed quasar, with a flux density ratio of 17.0 ± 0.6. 
C and D are probably lobes situated either side of the quasar 
core which lie too far from the lensing galaxy to form multiple 
images. The off-source rms is 0.09 mjy/beam and the beam is 
0'.'20 X O'.'IS at a position angle of 35°. The co-ordinates are offset 
from 03 16 50.933 +43 28 19.31, the greyscale is in mjy/bcam, 
and the contours increase in factors of 2 from 0.5 mjy/beam. 

The maps of J0316+4328 with the VLA and MERLIN 
are presented as Figures [1] [2] and [3] while fits to the compo- 
nents are presented in Table [2] We also mapped the system 
in Stokes Q, U and V, using the MERLIN data, and de- 
tected no polarised emission at an rms of 0.13 mjy/beam. 
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Table 2. Gaussian fits to the components of J0316+4328, giving positions relative to component A in the 22 GHz VLA map, integrated 
flux densities and sizes (deconvolved major axes). Components C and D were clearly resolved in the 6.0 GHz MERLIN map (see Figurc[3jl, 
so in this case 2 gaussians were fitted. The flux density is the total flux density of the 2 components, the size is their separation. B is 
unresolved in each map, while A is marginally resolved with MERLIN. The off-source rms values were 0.13, 0.09 and 0.09 mjy/beam 
and the beam sizes were 0'.'051 X 0'.'045, 0'.'20 X 0'.'18 and O'.'ll X 0'.'07 in the 6.0, 8.5 and 22.5 GHz maps, respectively. 
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Figure 2. 22.5 GHz VLA A array map of J0316+4328. Com- 
ponents A and image B are both compact, with a flux density 
ratio of 19.6 ± 1.3. With higher resolution, components C and 
D are clearly extended. The off-source rms is 0.09 mjy/beam 
and the beam is O'.'ll X 0'.'07 at a position angle of 25°. The co- 
ordinates are offset from 03 16 50.933 +43 28 19.31, the greyscale 
is in mjy/beam, and the contours increase in factors of 2 from 
0.3 mjy/beam. 



Figure 3. 6.0 GHz MERLIN map of J0316+4328. Component 
B is compact, component A is marginally extended (major axis 
0'.'051 at a position angle of 54°). The flux density ratio of compo- 
nents A and B is 19.2±0.8, taking the total flux density of A. Com- 
ponents C and D are resolved by MERLIN. The off-source rms is 
0.13 mjy/beam and the beam is 0''051 X 0''045 at a position angle 
of 61°. The co-ordinates are offset from 03 16 50.93 +43 28 19.3, 
the greyscale is in mjy/beam, and the contours increase in factors 
of 2 from 0.4 mjy/beam. 



J0316+4328 has 2 compact components (A and B) with flux 
density ratios of 19.6 + 1.3, 17.0 + 0.6, 19.2 + 0.8 at 6.0 GHz, 
8.5 GHz and 22 GHz, respectively. Errors are 1-a values 
taken from the map rms. The 2 steeper spectrum compo- 
nents (C and D) either side of A and B are extended with 
the VLA and clearly resolved by MERLIN. These are proba- 
bly lobes situated either side of the radio source's core which 
lie too far from the lensing galaxy to be multiply imaged. 

We observed J0316+4328 with the Very Long Baseline 
Array (VLBA) at 8.4 GHz on 2006 July 16. The observation 
was made at an observing rate of 256 Mb/s and included 113 
minutes on J0316+4328. Self-calibration was used to derive 
the phase solutions, and it was not possible to fit good solu- 
tions on the longest baselines, so the Mauna Kea, St. Croix 
and Hancock antennas were omitted. Components A and B 
(Figures |3] and O were found to have reduced flux densities 
compared to those in the VLA maps (Table [3}, presumably 
due to resolving out extended emission. No emission was de- 



tected at the locations of components C and D; they were 
completely resolved out by the VLBA. Component A has 
84% of its flux density in a central point source Al and the 
remainder in jet components A2 and A3. Component B has 
a single point source sub-component Bl, with a flux density 
ratio relative to A 1 of 23.7 ± 1.5, matching the 6 GHz and 
22 GHz ratios at the 2-<r level. 

Some simple lens models which fit the bright images 
place the central image only ~ 20 mas from image B. 
While these models are severely underconstrained (due to 
the dearth of information on the lens galaxy), it is plausi- 
ble that the VLA and MERLIN maps would blend image B 
with the central image, meaning that only the VLBA maps 
constrain the central image. No source was seen between 
images A and B, with a 5-cr limit of 0.33 mjy, although in- 
terpreting a non-detecti on over many resolu tion elements is 
not straightforward, see IZhang et all (|2007m . 
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Table 3. Gaussian fits to the components of J0316+4328 in the VLBA maps, giving positions relative to Al, integrated flux densities 
and sizes (deconvolved major axes). The blank field rms values were (0.08,0.065) mjy/beam in the A and B fields, respectively. The 
beam size was 2.8 X 2.5 mas. Al, A2 and Bl are compact. 
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Figure 4. 8.4 GHz VLBA map of the A component in 
J0316+4328, omitting the Mauna Kea, Hancock and St. Croix 
antennas. At VLBI scales component A splits into an unresolved 
core Al and jet components A2 and A3 at distances of 7 and 
43 mas and position angles of 48° and 53°, respectively. This 
matches the 51 mas size and 54° position angle of component A 
in the 6.0 GHz MERLIN map (Figure EJ. The off-source rms is 
0.08 mjy/beam and the beam is 2.8 X 2.5 mas at a position an- 
gle of 16°. The co-ordinates are offset from component Al, the 
greyscale is in mjy/beam, and the contours increase in factors of 
2 from 0.4 mjy/beam. 



3 DISCUSSION 

Although we lack high-resolution optical data, we are very 
confident that J0316+4328 is a gravitational lens. Compo- 
nents A and B are very likely to be lensed images of the same 
quasar, as they are compact at resolutions from 0'.'2 to 0'.'05 
and exhibit the same flux density ratio over a factor of 3.5 
in frequency (Figure [6]). At 3 mas resolution substructure 
is seen only in the brighter, more magnified image. All but 
one of the CLASS candidates with multiple co mpact compo- 
nents in VLBA maps were confirmed as lenses (jBrowne et al.l 
120031 ). The exception, B0827+525, had significantly differ- 
ent spe ctra for the 2 componen ts and is likely to be a binary 
quasar llKoopmans et al.ll200d ). 

We plan to obtain a VLBI map with rms ~ 0.02 mjy, 
which would definitively confirm the lensing interpretation 
by detecting parity-reversed substructure in image B match- 
ing that in image A. If B is a lensed image 19 times fainter 
than A, it should include sub-components B2 and B3 with 
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Figure 5. 8.4 GHz VLBA map of the B component in 
J0316+4328, omitting the Mauna Kea, Hancock and St. Croix 
antennas. It appears as a single compact point source Bl, with a 
flux density ratio relative to Al of 23.7 ± 1.5. The off-source rms 
is 0.065 mjy/beam and the beam is 2.8 X 2.5 mas at a position 
angle of 16° . The co-ordinates are offset from the brightest com- 
ponent Al (shown in Figure |4j, the greyscale is in mjy/beam, 
and the contours increase in factors of 2 from 0.2 mjy/beam. 



flux densities of 0.08 and 0.16 mjy, and separations from Bl 
of 1.7 and 10 mas. This map would also improve the flux 
density measurement of Bl (which currently limits the pre- 
cision of the VLBI flux density ratio) , and might even detect 
a central image. 

We are also pursuing high resolution optical images 
with the William Herschel telescope and spectroscopic ob- 
servations with the Lick telescope. The USNO Bl catalog 
shows an optical object with R=19.1 and B=20.3 at 03 16 
50.8693 +43 28 19.880, an offset of -Of! 69 in right ascension 
and +0'.'55 in declination from image A in the MERLIN 
map. The USNO c atalog has astrometric accuracy of 0'.'20 
|Monet et al.l [2003) and this object has errors of / .'38 in 
right ascension and O'.'lO in declination, while the MERLIN 
map has astrometric accuracy and positional error both bet- 
ter than 0'.'05. The optical and radio positions are consistent 
at the 2-<t level, but it is not possible to identify the USNO 
object as a lensed image, the lens galaxy, or a blend of both. 
High resolution optical images should confirm the lensing 
interpretation, and allow detailed lens modelling. 

For isothermal profiles with moderate ellipticity or 



Asymmetric Double Lens J0316+4328 5 



E 





1 r 


1 1 — 


i 

A 


■ 1 






• 


C 




















D ' 




























• 












i 





10 

Frequency (GHz) 



Figure 6. Spectra of the components of J0316+4328 from the 
VLA and MERLIN observations. A and B (solid squares, solid 
lines) have the same spectrum from 6 GHz to 22 GHz, while C 
and D (hollow squares, dashed lines) are noticeably steeper. The 
VLBA flux densities of Al and Bl are plotted as solid circles. 
Component A and Al flux densities are multiplied by 0.25. 



4 CONCLUSION 

We have detected a probable gravitational lens with 2 bright 
images at a considerably higher flux density ratio than any 
previously known double lens. Improved VLBI and high- 
resolution optical observations are underway. This highly 
asymmetric lens should give unusually stringent constraints 
on the density profile of the lens galaxy. 
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shear, the lensing cross-section for asymmetric doubles 
(magnification ratio 10 — 20) is ~ 20% that of quads and 
symmetric doubles (magnification ratio 1 — 10) combined. 
With 21 quad or symmetric double lenses in CLASS, 4-5 
asymmetric doubles would be expected. 3 such systems have 
been found, counting the confirmed lens B1030+074 and the 
2 strong candidates J0316+4328 and J0935+0719. 

Our lens modelling must accommodate the extended, 
steep-spectrum components C and D. Their radio spectra 
are different from each other as well as from A and B (Fig- 
ure |6}, and their separation is 1'.'12 as opposed to 0'.'40 for 
images A and B, so C and D are unlikely to be a pair 
of lens images. Their orientation either side of image A, 
and the fact that substructure in A points towards the 
brighter component C, argues that they are diffuse lobes 
lying either side of the source's core, in a typical double ra- 
dio source configuration. With the lens image separation of 
C//40 arguing for a relatively less massive lens galaxy with 
a small lensing cross-section, it should be possible to find 
models where the source's lobes lie outside the multiply- 
imaged region of the source plane, while the core lies just 
within it. The configuration actually resembles that of the 
first lens system B0957+561 (IWalsh. Carswell. fe Wevmannl 
1 19791 ; iRoberts. Greenfield, fc Burkdll979h . 

Our modelling will give strong constraints on the lens 
galaxy's density profile. For near isothermal models which 
are a good approximation to most lens galaxies (|Rusin et al.l 
l2003l ; lKoopmans et al.ll2006h . images A and B form at radii 
differing by a factor of 15-20 and constrain the density pro- 
file over this large range of galactocentric radius. The con- 
straints will be particularly strong if substructur e is detected 
in ea ch image, as in the case of B1152+199 (|Rusin et al.l 
2002). Also, central images become brighter as image B be- 
comes fainter, relative to image A. Compared to quad or 
symmetric double lenses, J0316+4328 is more likely to show 
a central image in sensitive VLBI maps, and even an up- 
per limit on the central image flux density will be a more 
stringent model constraint. 
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